Molecular prototypes for spin-based CNOT quantum gates 
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We show that a chemically engineered structural asymmetry in [Tb2] molecular clusters renders 
the two weakly coupled Tb 3+ spin qubits magnetically inequivalent. The magnetic energy level 
spectrum of these molecules meets then all conditions needed to realize a universal CNOT quantum 
gate. 

PACS numbers: 75.50.Xx,03.67.Lx,75.40.Gb 



Quantum computation^ 2 - relies on the physical realiza- 
tion of quantum bits and quantum gates. The former can 
be in any of two distinguishable and addressable states, 
denoted here as spin up | ff) and spin down | J|), but 
also, as opposed to their classical counterparts, in any 
arbitrary linear superposition of these. The latter in- 
volve controlled operations on two coupled qubits X The 
universal controlled-NOT (CNOT) gate is the archetype 
of such a controlled operation. It flips the target qubit 
depending on the state of the other, or control qubit (see 
Fig. QJa)). In terms of its material realization, this defi- 
nition implies that each of the two qubits should respond 
inequivalently to some external stimulus, e.g. electric or 
magnetic fields. 

Solid-state candidates to build these basic ingredi- 
ents include superconducting circuits^"— single spins 
in semiconductors and, more recently, molecular 
nanomagnets^ — The latter are attractive for scalability, 
since macroscopic arrays of identical molecular qubits can 
be prepared via relatively simple chemical methods i 15 i 16 
State-of-the-art achievements include the measurement 
and minimization of decoherence ratesii — and the syn- 
thesis of mutually interacting qubit pairs j 2Q i 21 However, 
the realization of a universal CNOT quantum gate inside 
a molecular cluster remains an outstanding challenge. 14 
Here, we show that asymmetric [Tb2] molecular clusters 
show also the magnetic asymmetry that is required for 
the realization of spin-based CNOT gates. 

Lanthanide ions are promising candidates for encoding 
quantum information, 22 The magnetic energy levels are 
split by the crystal field generated by surrounding atoms, 
which often results in a well-isolated ground state dou- 
blet. Each metal ion behaves then as an effective (two- 
state) spin- 1/2, therefore providing a good realization of 
qubit states. For the realization of a spin-based qugate, 
it seems therefore natural to look for simple molecules, 
made of just two weakly-coupled lanthanide qubits. How- 
ever, the synthesis of asymmetric molecular dimers is not 
straightforward, as Nature tends to make them symmet- 
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FIG. 1: (Color online) (a) Schematic illustration of a quantum 
CNOT operation on two coupled spin qubits. (b) Targeted 
synthesis^ of an asymmetric [Tb2] complex using three H3L 
ligands (depicted in green or black colors in the complex), 
(c) The resulting N2O6CI and N3O6 coordination polyhedra 
around the two Tb 3+ ions strongly differ, exhibiting Ci sym- 
metry and C^v symmetry, respectively (see [24| for details), 
and induce a misalignment between the two local anisotropy 
axes. Color code: "control" Tb, purple; "target" Tb, dark 
blue; CI, orange; O, red; N, light blue. 
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ric. We propose a simple solution that exploits the ability 
of chemical design to finely tune the internal molecular 
structure. The chemical reaction is sketched in Fig.QJb). 
By the use of an adequate asymmetric organic ligand^ 
H3L, we synthesize a dinuclear complex of Tb 3+ ions, 
hereafter briefly referred to as [Tb2], in which the metal- 
lic dimer is wrapped by three ligands. Therefore, each 
metal ion in the molecule is in a different coordination 
environment (see Fig. []Jc)). In the following, we de- 
scribe physical experiments that enable us to verify that 
this molecular cluster fulfills the basic conditions to act 
as a CNOT, namely the appropriate definition of the two 
qubits, the existence of a weak coupling between them, 
and the magnetic asymmetry. 

Magnetic and calorimetric measurements were per- 
formed on powdered specimens. Above 1.8 K, ac (h ac = 
4 Oe) susceptibility and magnetization data were mea- 
sured with a commercial SQUID magnetometer. Ac 
susceptibility data were also measured with a /iSQUID 
susceptometer^- operating in the 13 mK < T < 1.5 K 
temperature and 0.01 Hz up to 1 MHz frequency ranges. 
The amplitude of the ac magnetic field was, in this case, 
h ac = 1 mOe. Dc-magnetization measurements below 
2 K were performed using a homemade Hall microprobe 
installed in a dilution refrigerator. Heat capacity mea- 
surements down to ~ 0.35 K were performed using the 
relaxation method by means of a commercial setup. 

The ac magnetic susceptibility x provides direct in- 
sight on the magnetic anisotropy of the Tb 3+ ions and 
their mutual coupling. For the lowest frequency (0.0158 
Hz) and above 100 mK, the in-phase component x' 
gives the equilibrium paramagnetic response (see Fig. 
[2] (a)). In this regime, the cluster's effective mag- 
netic moment /i e ff can be approximately determined as 
Meff — (3/cbX^/^a) 1 / 2 - At room temperature, /i e ff = 
13.7(l)/iB agrees with the effective moment of two uncou- 
pled free ions, i.e. /i e ff = #jMb[2J(J + l)] 1 ^ 2 = 13.74/iB, 
where gj = 3/2 and J = 6 are the gyromagnetic ra- 
tio and the total angular momentum given by Hund's 
rules. The drop observed below approximately 100 K 
can be assigned to the thermal depopulation of mag- 
netic energy levels split by the crystal field. The value 
/i e ff = 12.5(l)/iB measured between 3 and 10 K is close to 
the "Ising" limit /i e ff = #jMb2 1 / 2 J = 12.72/iB, character- 
istic of two uncoupled Tb 3+ ions whose angular momenta 
Ji and J2 can only point either up or down along their 
local anisotropy axes. The magnetic anisotropy can be 
determined by fitting in the intermediate tempera- 
ture regime, using expressions derive d 24 ! 26 for the sim- 
plest uniaxial anisotropy T~L = —D(jf z + J$ z )- The fit 
gives D/kB = 17 K. For such a strong anisotropy, excited 
levels are separated by more than 180 K from the ground 
state doublet, associated with the maximum projections 
mj = zb J (see Fig. EJb)). These states provide then, for 
each ion, a proper definition of the qubit basis | ft) and 
I Jj), namely: 




FIG. 2: (Color online) (a) Ac susceptibility of [Tb2] measured 
on a powdered sample for different frequencies of the alternat- 
ing excitation field. Main panel: x'T product (left axis) from 
which /i e ff is determined (right axis). Inset: x (solid symbols) 
and x" (open symbols). Dc susceptibility data measured un- 
der a constant magnetic field /noH = 0.1 T are also shown. 
The lines represent least square fits of the ac (solid lines) and 
dc (dotted lines) susceptibilities based on the spin Hamilto- 
nian (|3]) for collinear (6 = 0, red thin lines) and noncollinear 
(5 = 66°, blue thick lines) anisotropy axes, (b) Zero- field 
energy level structure of [Tb2] derived from these fits. 



\t) = \J = 6,m J = 6) (1) 
\^) = \J = Q,mj = -6) (2) 

By further decreasing temperature below 3 K, we ob- 
serve a second drop in x'T (Fig. [2j) that we associate 

with the antiferromagnetic coupling — 2J ex JiJ2 of the 
two qubits. This interpretation is corroborated by the 
magnetic heat capacity c m> p, shown in Fig. [3] (see I24]). 
At zero field, c m ,p shows a Schottky-type broad anomaly 
centered at T max ~ 0.9 K, which arises from the en- 
ergy splitting e between antiferromagnetic (| ft) ® | 4) 
and I 4) ® I ft)) and ferromagnetic (| ft) <8) | ft) and 
I J|) I 4)) states. Using the condition /cBT max = 0.42e, 
we determine the splitting e/k^ — 4:J ex J 2 /kB — 2.14 K 
and, from this, the effective exchange constant J ex /^B = 
-0.016(1) K. 

Considering the value of e, the monotonic increase of x' 
down to 100 mK would be puzzling unless the magnetic 
moments of the two Tb 3+ ions do not exactly compen- 
sate each other. As the coordination sphere determines 
the magnetic anisotropy, the anisotropy axes of the two 
ions need not be parallel to each other, but can instead 
make a tilting angle J, as shown in Fig. [TJ Because of 
this misalignment and the very strong anisotropy, even 
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FIG. 3: (Color online) Magnetic heat capacity of [TD2]. (a) 
Data measured on a powdered sample under different values of 
the applied magnetic field, (b) and (c) Theoretical predictions 
for collinear and noncollinear anisotropy axes, respectively. 



the states (| ff) (g) | J|) and | J|) <g> | ff)) preserve a net 
magnetic moment. More importantly, the two ions will 
couple differently to an external magnetic field, i.e. their 
effective gyromagnetic ratios g\ and #2 will be different. 
For instance, if H is applied along one of the anisotropy 
axes, say of qubit "1", g\ = gj whereas #2 = #jcos£, 
i.e. the field makes the two spins inequivalent, just as 
required for the CNOT operation. 

Susceptibility and heat capacity measurements confirm 
that gi 7^ #2 • The absence of magnetic cancelation in the 
cluster causes the paramagnetic response observed be- 
tween 10 K and 100 mK (see the inset of Fig. [2]). Indeed, 
below 10 K x'T (and thus also /i e ff) is much larger than 
predicted for collinear anisotropy axes (5 = 0). In con- 
trast, an excellent agreement is obtained for 5 = 66°. We 
have also measured the dc-susceptibility under a mag- 
netic field of 0.1 T, obtaining further experimental evi- 
dence for noncollinear axes with the same values of J ex 
and 5 estimated above. 

The same conclusion is drawn from heat capacity data 
measured under nonzero applied magnetic fields, shown 
in Fig. 03 This quantity reflects the magnetic field de- 
pendence of the energy level structure, which, in its turn, 
should strongly depend on S. As with the magnetic data, 
the results are in qualitative and quantitative agreement 
with the calculations made for 5 = 66°. 

As an independent confirmation of the magnetic in- 
equivalence of the two Tb 3+ ions, we have also measured 
magnetization isotherms at fixed temperatures T = 0.26 



and 2 K. The results are shown in Fig.|4^a). The presence 
of a non-compensated /i e ff has dramatic consequences 
at very low temperatures and for low magnetic fields, 
when the thermal populations of the excited ferromag- 
netic states (| ff) (g) I ff) and | J|) (g) | become negligibly 
small. For 5 = 0, the ground state of the cluster would 
be nonmagnetic, therefore M c± until it would abruptly 
change at the magnetic field hqH ~ 0.15 T where the 
energy level crossing occurs. Experimentally, we find 
instead a finite paramagnetic response starting already 
from H = 0, which shows that the molecular ground 
state possess a net magnetic moment (and therefore that 
S must be different from and close indeed to 66°). The 
level crossing does still give rise to an abrupt magneti- 
zation jump at a somewhat higher field, confirming the 
antiferromagnetic character of the exchange interactions 
and the value of J ex determined from heat capacity ex- 
periments. 
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FIG. 4: (Color online) (a) Magnetization isotherms of [Tb2] 
measured at two temperatures. The data at T = 0.26 K 
are compared with calculations made for collinear (6 = 0, 
dotted line) and noncollinear (5 = 66°, solid line) anisotropy 
axes, (b) Field-dependent magnetic energy levels of [Tb2] 
calculated by numerical diagonalization of the Hamiltonian 
for S = 66° . Only levels corresponding to magnetic nuclear 
states m/,i — raj, 2 = —3/2 are shown. 

All experiments presented so far lead us to conclude 
that the essential physics is captured by a Hamiltonian 
containing the uniaxial anisotropy, the exchange cou- 



plings, the Zeeman energy and, finally, the hyperfme in- 
teractions with the nuclear spins I = 3/2 of Tb. In addi- 
tion, we have already anticipated that the relevant qubit 
states | ff) and | J|) are related to mj = ±6. In this 
reduced subspace, the Hamiltonian simplifies tc 24 

T~L = — 2J e *Ji,zJ2,z — QiI^bHJi^ — g2^BHJ2,z 

+ Aj (J liZ I liZ + J 2 ,zl2,z) (3) 

where Aj /&b = 2.5 x 10 -2 K is the hyperfine constant. 

The Hamiltonian (j3j) enables us to finally discuss the 
performance of [Tb2] as a CNOT qugate. The ensuing 
energy level spectrum is shown in Fig. ffl^b) for 5 = 66°. 
For the sake of clarity, we show only levels associated with 
a single nuclear magnetic state of the cluster, namely 
with nuclear spin projections m/ 5 i = mi^ = —3/2. The 
magnetic asymmetry induced by the tilting angle S = 
66°, enables to univocally single out any of the desired 
transitions. For instance, at fi^H = 0.28 T, only the two 
states | ff)i (g) | ft) 2 and | ft)i (S> | 2 are separated by 
an energy equivalent to that of v = 9.8 GHz photons 
(X-band EPR). The existence of multiple nuclear spin 
states does not invalidate the proposed operation scheme 
because of the Am/ = selection rule. 

Concluding, we have shown that molecular clusters 
containing two lanthanide ions meet the ingredients re- 
quired to implement a CNOT quantum gate. The defini- 
tion of control and target qubits is based on the magnetic 
inequivalence of the two ions, which has been achieved by 
chemically engineering dissimilar coordination spheres. 



Although we have restricted here to discussing the [Tb2] 
compound, for which the magnetic asymmetry can be 
more easily determined on account of its large angular 
momentum, the same molecular structure can be real- 
ized with other members of the lanthanide series^ This 
flexibility, characteristic of molecular materials, enables 
a vast choice of qugate designs. For instance, it will be 
possible to select isotopes with zero nuclear spin (like 
most of Er isotopes) or ions with dominantly easy-plane 
anisotropy that, like Ce (J = 5/2), have a mj = ±1/2 
electronic ground state doublet. In the latter case, the 
CNOT operation would correspond to an allowed electro- 
magnetic transition, which should enhance the attainable 
Rabi frequencies to values of order 120 MHz/mT. Chem- 
ically engineered molecular qugates can therefore open 
promising avenues for the realization of scalable quan- 
tum computing architectures. 
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